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Abstract-An ideal multi-objective optimization method for economic emission 
load dispatch (EELD) with non-linear fuel cost and emission level functions in 
power system operation is presented. In this paper, the problem treats economy, 
emission, and transmission line security as vital objectives. The load constraints 
and operating constraints are taken into account. Assuming goals for individual 
objective functions, the multi-objective problem is converted into a unique- 
objective optimization by the goal-attainment method, which is then taken care 
of by the simulated annealing (SA) technique. The solution can offer a best 
compromising solution in a sense close to the requirements of the system 
designer. Results for a 30-bus IEEE test case system have been utilized to 


demonstrate the applicability and authenticity of the proposed method. 


I. INTRODUCTION 


Looking at the sophistication of the power utility 
sectors, basically for thermal power plants, the security 
measures for the power system as a whole are taken into 
account incorporating environmental effects out of 
generating units with economic aspects. The basic 
objective of economic emission and load dispatch[1-2] is 
to trace out the optimal power generated in fossil-based 
generating units by optimizing the fuel cost attaching a 
squared nonlinear dependence in the cost function and a 
non-smooth emission level function simultaneously 
taking into account various inequality constraints. Zahavi 
and Eisenberg [3] presented a method to solve economic 
environmental power dispatch problem without 
exemplifying it. Similarly In goal programming [4], 
approach for EELD, the transmission loss and the line 
security measures were not taken into account. 

A better multi-objective optimization procedure 
based on probability security criteria to obtain a set of 
non-inferior solutions was presented in [5]. Nanda et al. 
[6] have formulated the EELD problem with line flow 
constraints and solved it through a classical technique, 
but the mathematical formulation of the security 
constrained problem would require a very large number 
of constraints to be considered. This classical technique 
introduced a preference index between the two objectives 
(economy and emission) in order to decide on an optimal 
solution, and this would result in complex problem 
formulation when the number of objective functions 
exceeds two. The major disadvantage of the aforesaid 
methods in solving the EELD problem is that they are 
insufficient for handling non-smooth fuel cost and 
emission level functions. In other words, they use 
approximations to restrict severity of the problem. The 
insufficient accuracy induced by these approximations 1s 
not desirable. Simulated annealing (SA) can ameliorate 
this undesirable characteristic by simulating the physical 
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annealing process for the computation of the global or 
near-global optimum solutions for optimization 
problems. Amongst other applications, the SAtechnique 
has been successfully applied to economic dispatch [7] 
and hydrothermal scheduling [8]. 

As described in the synopsis above the EELD problem 
comprising the aforesaid objectives is converted into a 
single objective optimization problem using goal 
attainment (GA) method which is later on dealt by 
Simulated Annealing (SA) technique for seeking 
reasonably approximate global optimum solution in spite of 
presence of non-smooth unit characteristics. The algorithm 
developed has been implemented on 30-bus power system 
consisting of 5 thermal generators and 41 transmission 
lines. The experimental results are also presented. 


Il. FORMULATION OF THE PROBLEM 


The present formulation uses the Economic 
Emission Load Dispatch problem as a multi-objective 
optimization problem which is concerned with an attempt 
to optimize each objective simultaneously. Care is taken to 
see that the equality and inequality constraints of the 
system are satisfied. The following objectives and 
constraints are taken into consideration in the formulation 
of the Economic Emission Load Dispatch problem. 


A. Objectives 

Economy. Consider a system having N buses and 
NL lines. Let the first NG buses have sources for power 
generation. Taking into account the valve-point effects[9], 
the fuel cost function of each generating unit is expressed 
as the sum of a quadratic and a squared sinusoidal function. 
Therefore, the total cost of generation C in terms of control 
variables PG ’s is given by the following expression: 


NG 
f,(PG) = C =>°0.5a,PG? +b,PG, +c, +|d, xsin’(eix(PG" — PG,))|$/h (1) 
i=1 


where PG; is the real power output of an ith generator, NG 
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is the number of generators, and a;, b;, c;, d;, e; 18 fuel cost 
curve coefficients of an ith generator. 


Emission. The power generating stations being the primary 
sources of nitrous oxides, they are strongly objected by the 
Environmental Protection Agency to reduce _ their 
emissions. In this study, nitrous oxide (NO,) emission is 


NG 
f,(PG)=E=)'0.5a,PG; + BPG, +, +, exp(k,PG,) Ib/h (2) 


i=] 


where @,, {,,7,,77, and k, are emission curve coefficients 


of the ith generating unit. 


B. Line Security. 

Security constraints involve critical lines for 
replacing huge no of transmission lines in the power 
system network that are of immense importance in deciding 
the optimal solutions for an electric power systems. The 
system designer interprets the transmission lines violating 
the equality and inequality constraints as critical lines. The 
security constraints of the system can give better prospects 
by optimizing the following objective function: 


fy(PG) =S = DL j(PG)/E™) (3) 


where L,(PG)is the real power flow, Lis the 


maximum limit of the real power flow of the jth line, and k 
is the number of monitored lines (critical lines). The line 
flow of the jth line is expressed in terms of the control 
variables PG ’s by utilizing the Generalized Generation 
Distribution Factors(GGDF) [11], and is given below: 


F (PG) = s. (D ,PG,) (4) 


i=l 


where D j_ 18 the generalized generation distribution factor 


(GGDF) for line j due to generator i. 


Load Constraint. The real power balance between 
generation and the load is maintained always thinking the 
load at any time to be constant: 


NG 

> PG, =P, +P, (5) 

i=l 
where P,, is the total real power demand and P, is the 
total real power loss. The latter is represented as [10] 


P= > APG | (6) 


where A, is the loss coefficient due to the generator i. 


The loss coefficient are evaluated from base load flow 
solution. 


Operating Constraints. For achieving stable operation each 
generating unit is to be confined within its lower and upper 
real power limits. 

PG™ <PG, < PG™ (7) 


where PG. and PG," are the minimum and 


maximum real power output of ith unit, respectively. 
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taken as the selected index from the viewpoint of 
environment conservation. The amount of emission from 
each generator is given as a function of its output [8], 
which is the sum of a quadratic and an exponential function 
in the present work. Therefore, the total emission level E 
from all the units in the system can be expressed as 


Hl. THE GOAL-ATTAINMENT METHOD 


Multi-objective formulation index is dealt with a 
set of objectives f(x) = [fi(x), fa), ......... , Jn (x)].In this 
method the designer sets a vector of designed goals g 
=[21,2, «2, |? which form a powerful tool[13-15]that 
associate with aforesaid objectives.The level of attainment 
of the goals is controlled by a weight vector w = [w,wo, 
siay WH | 
In this GA method of optimization , the aforesaid nonlinear 
problem is solved: 


(8) where 


x is a set of desired parameters which can be varied, A is 


minimize 2 x €Q subject to g+ Aw = f(x), @EA,. 


a scalar variable which introduces an element of slackness 
in to the system, W is a feasible-solution region that 
satisfies all the parametric constraints, and 


A. ={@E R'St. @, 26> a, =landez0 } 


i=] 


Figure | illustrates two-dimensional goal-attainment 
method. The multi-objective optimization is concerned 
with the generation and selection of non-inferior solution 
points [15] to characterize the objectives where an 
improvement in one objective requires a degradation in the 
others. 





w g 
1 1 


(x) 
Le , 


Figure 1: Illustration of two-dimensional goal-attainment 
method. 


By varying w over A. , the set of non-inferior 


solutions is generated. In the two _ dimensional 
representation of Figure 1, the set of non-inferior solutions 
lies on the curve AB. The weight vector w enables the 
designer to express a measure of the relative trade-offs 
between the objectives. Given the vectors w and g, the 


direction of the vector g + Aw is determined. A feasible 
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point on this vector in function space which is closest to the 


origin is then searched. The first point Ay at which g + 


A w intersects the feasible region F in the function space 
would be the optimal non-inferior solution. During the 


optimization, ZA is varied, which changes the size of the 
feasible 

region. The constraint boundaries converge to the unique 
solution point 


{fio(x), foo(x)}.For optimizing the value of A the Simulated 
Annealing method is employed as under: 


IV. SIMULATED ANNEALING TECHNIQUE 


In this method[16-17] a candidate solution is 
generated which is accepted when it becomes a better 
solution to generate another candidate solution. If it is 
deteriorated solution, the solution will be accepted 
provided its probability of acceptance P,(A) given by 
equation (9) is greater than an arbitrarily generated number 
between 0 and 1. 

P(A) = [1/{1 + exp(A/T)}] (9) 

where A is the amount of deterioration between the new 
and the current solutions and 7 is the temperature at which 
the new solution is generated. In forming the new solution 
the current solution is perturbed [5] according to the 
Gaussian probability distribution function (GPDF). The 
mean of the GPDF is taken to be the current solution, and 
its standard deviation is given by the product of the 
temperature and a scaling factor+. The value of + is less 
than one, and together with the value of the temperature, it 
governs the size of the neighborhood space of the current 
solution and hence the amount of perturbation. The new 
solution is formed by adding the amount of perturbation to 
the current solution.In the next iteration the temperature is 
reduced according to a cooling schedule. The following 
geometric cooling schedule is adopted in the present work 
[15]: 

Ty= 2 P To (10) 

where 7) and 7), are the initial temperature and the 
temperature at the vth iteration, respectively, and r is the 
temperature reduction factor. The solution process 
continues until the maximum number of iterations is 
reached and the optimum solution is found. 


V. EXPERIMENTAL RESULTS 


The algorithm developed in the previous section 
has been applied to a 30-bus test system. The system 
consists of 5 generators and 41 lines. The line data and the 
load data are given in the Appendix. Table 1 gives the real 
power operating limits whereas Tables 2 and 3 give the 
cost curve and emission curve coefficients of the 5 
generators. The voltage at the 5 buses are kept fixed 
respectively, to the values 1.0634,1.0482,1.0354,1.008 and 
1.0631 p.u. .The system load was taken on IOOMVA base. 


In applying the developed algorithm for the test 
system, the appropriate values of the control parameters are 
set. These parameters are initial temperature 7, the scaling 
factor + for GPDF, the temperature reduction factor r, 
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Figure 2: FLOW CHART FOR Multi-objective Economic 
Emission Load Dispatch with Nonlinear Fuel Cost and 
non-inferior Emission Level Functions 


maximum number of iterations VMAX, and the number of 

trials per iteration TMAX. In the present work 7p, +, 
VMAX, and TMAX were set, respectively, to the values of 
50000, 0.02, 200, and 2000. As per a guideline [17], the 
value of r lies in the range from 0.80 to 0.99. SA with a 
slow cooling schedule usually has larger capacity to and 
the optimal solution than with a fast cooling schedule. 
Hence, for seeking the optimal solution the value of r is 
required to set close to 0.99 so that a slow cooling process 
is simulated. The appropriate setting of r was set by 
experimenting its value in the range from 0.95 to 0.99, and 
this value was found to be 0.98. The following different 
case studies were conducted to illustrate the performance of 
the proposed algorithm. 


The variation of cost and emission functions with 
the optimized real power generation for the five test case 
systems is illustrated in plot-1 and plot-2 respectively. 
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Figure 3: Variation of cost function 


f, (PG, ) with PG, (plot-1) 
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Figure 4: Variation of emission level function f, (PG, ) 


with PG, (plot-2) 


The variation of cost and emission functions with the 
optimized real power generation for the seven test case 
systems is illustrated in plot-3 and plot-4 respectively. 
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Figure 5: Variation of cost function f,(PG,) with PG, 
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Figure 6: Variation of emission level function f, (PG; ) 


with PG, (plot-4) 
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This shows a better computation result as the 
number of buses get extended at the cost of little high time 
of computation. However the computation time can be 
reduced further by using advanced computing machine and 
differential evolutionary technique algorithm for multi- 
objective minimization. 

Table 1: Operating limits (p.u.) of generators on 
100 MVA base 


Gen | PGmin; | PGmax; 


[sa 


Table 2: Cost curve coefficients of generators 
Ci d; ej 


[ns 00 [ao [30 [om 
2 fo.0030 | 7.8000 | 60.0 [140 |0.040- 
3 [0.0012 | 2.1000 | 100.0 | 160 | 0.038 


Emission curve coefficients of generators 








Table 3: Cost curve coefficients of generators 


ete patriye 
# 1 
2 | 010030 | 7.8000 | 60.0 | 140 [0.040 





Case 1 

A goal vector g = [g1,82,g93]’ = [357.1, 228.05, 
1.1503]was generated automatically using Step 1 (iv) of the 
computational algorithm. Here, g; is the generating cost 
objective being expressed in $/ h, g. is the emission level 
objective being expressed in Ib/ h, and g31s the line security 
objective. The vector w = [wy,w2,w3]’= [0.3, 0.5, 0.5]’ that 
signifies the preference direction of the designer towards 
the goals was given by the designer. 
Case 2 

A goal vector of g = [389, 251.06, 1.61]? along with an 
identical weight vector as in case 1 is considered in this 
case. The goal values were assumed given by the Designer 
through his/ her experiences. 

The optimum generation schedule was determined for 
both cases and are presented in Table 4. 

The following three additional cases are also 
considered in the present study to observe the closeness of 
the value of a particular objective function towards its goal 
Case 3 
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A goal vector of g = [399, 252, 1.90]’ along with a 
weight vector of w = [0.4, 0.3, 0.5]’ is considered. 
Case 4 

In this case, a goal vector of g = [400, 257, 2.100]’ 
and a weight vector of w = 
[0.5, 0.4, 0.5] are considered. 
Case 5 

This case considers a goal vector of g = [397, 227, 
2.30]’ and a weight vector of w = [0.5, 0.5, 0.4]’. 

Case 3, case 4, and case 5 were computed using the 
proposed algorithm with same control parameters as in the 
previous two cases. The optimum generation schedule and 
the values of objective functions for these three cases are 
presented in Table 4. 


The algorithm has been implemented in the above method 
using MATLAB programming language and the software 
system are run on a 2.53GHz computers. 


VI. CONCLUSION 


The simulated annealing method along with goal 
attainment method was used to solve the aforesaid EELD 
problem with non-linear cost and emission functions 
characteristics in function space. Specifically the squared 
value of sin term attached to the cost function minimize the 
generation cost as depicted by the result analysis through 
table 1-5. An advantage of the proposed method is that it 
does not impose any convexity restrictions on the 
generating unit characteristics. In addition, it also allows 
the Designer to decide on different preferences for the 
objectives toward the goals according to the system 
operating conditions, thus resulting in a more flexible 
operation on generating units. The only demerit of the 
proposed method is longer execution time that can be 
improved by further developing the algorithms and using 
advance processors for computation purpose. 
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